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TECHNICAL NOTE 4090 . 

ANALYSIS OF SHOCK MOTION IN DUCTS DURING 
DISTURBANCES IN DOWNSTREAM PRESSURE 
By Herbert G. Hurrell 


SUMMARY 

The effect of small downstream pressure disturbances on the position 
of a normal shock in a duct with area variation is analyzed. For t he 
analysis, the gas flow is treated as quasi- one-dimensional, nmri boundary 
layer is neglected. The analysis shows that there is a first-order lag 
relation between shock position and small downstream disturbances in 
pressure which occur at frequencies below a given limit. The time con- 
stant and the gain of this lag are expressed in terms of a dimensionless 
time constant that depends -only on the steady-state Mach number of the 
shock. 


INTRODUCTION 

When studying the dynamic behavior of propulsion systems for super- 
sonic flight, numerous problems arise that involve the motion of shock 
waves in ducts. In seme of these problems, the shock defines one bound- 
ary condition of the transient flow being analyzed. This condition 
occurs, for example , when the problem is one of determining the dynamics 
of inlet diffusers for control purposes. In other problems, it is the 
shock transient itself that is of primary concern. Such is the case when 
an inlet diffuser with severe buzz characteristics is involved. 

Several presentations of the basic theory of shock motion are to be 
found in the literature; reference 1 is, perhaps, the most comprehensive. 
In reference 2, however, the discussion is centered on the problems of 
duct flows. In this reference, a linearized equation governing shock 
motion is formulated for the case of a normal shock set off from its 
equilibrium position in an otherwise steady flow. The relaxation time 
for the return of the shock to equilibrium is thus determined and is used 
in a discussion of shock- wave stability in diverging or converging ducts. 

In the present report, a linearized analysis is made to define the 
transient imposed on a normal shock by an unsteady downstream flow. The 
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forcing variable of interest is the pressure downstream of the shock. 

The purpose is to present the dynamic relation between the "shock position 
and downstream, pressure in a form readily usable in the study of engine 
dynamics . 

In the analysis, which is made for small perturbations of the shock, 
the gas flow is treated as quasi- one -dimensional and boundary layer is 
neglected. The treatment of the unsteady flow downstream of the moving 
shock Is simplified by imposing a frequency "limit on the pressure disturb- 
ance. The relation sought is derived as a Laplace transfer function, 
which is considered the most convenient form for most uses . 

SYMBOLS 

A cross-sectional area of duct I 

a speed of sound 

c specific heat at constant pressure 

specific heat at constant volume 

D 3d 

substantial derivative, + u ^ 

k gain 

M Mach number 

P total pressure 

p static pressure 

S entropy 

s complex operator 

t time 

u gas velocity 

v shock velocity relative to duct 

x position coordinate (positive in downstream direction) 

X ratio of specific heats 
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p density 

ct time constant 

o' dimens ionless t ime constant 

cd frequency 

Subscripts: 

ss steady state 

r reference (or datum) values 

t total conditions (i.e., stagnation conditions if gas is brought 

to rest isentrqpically) 

0 entering shock 

1 leaving shock 

2 at fixed station just downstream of steady-state shock position 
Superscript: 

* critical conditions (i.e., conditions where gas velocity is equal 

to sonic velocity) 


AMLYSIS 

Because of its extreme thinness, a shock wave may be treated as a 
discontinuity with accuracy. With such a treatment, of course, the un- 
steady terms vanish from the fundamental relations across the shock. 

Flow conditions downstream of a moving shock, therefore, can be related 
to those upstream by steady-state equations, provided that the velocities 
used are relative to the shock. Hence, for a normal shock, moving at a 
velocity v with respect to the duct, the pressures on the two sides of 
it are related as follows: 



where aJ 3 velocities are considered positive in the downstream direction. 
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Consideration is new limited to small perturbations from the steady- 
state condition. When terms of higher order than these perturbations are 
neglected, equation (l) yields 


- rh: («o - + 


where only the A quantities, which represent the small perturbations, 
are time dependent; the coefficients of these variables are steady-state 
values. The amission of the higher order terms is tantamount to imposing 
the following conditions: 


Av 

Vo 

Mo 

Apo 

Pq 


« 1 


« 1 


« 1 


The gas flow upstream of the shock, of course, is not affected by a 
downstream disturbance . The variables AMq and Ap Q in equation (2), 

therefore, depend only on the area of the duct; and if the area variation 
is fixed, they are related to the shock position xq as follows: 


The velocity 
tion since 



of the shock can also be expressed in terms of shock posi- 


Av 


d(Ax 0 ) 

dt 


where the position coordinate x is positive for the downstream 
direction. 
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With the use of these expressions for APq, and Av, equation 

(2) becomes 


^>1 “ 


4r 

r + 1 




~(^o\ (M o' 
V^oJggV^o 


^0 - f- 

a o 


d(Ax 0 ) 


dt 


*1 


Axoy ^0 

ss' u/ ss 


(3) 

The shock movement Axq, therefore, is defined in terms of the pressure 
change Ap-^. This pressure variation occurs at the downstream face of 
the moving shock. The procedure now is to relate Ap-j_ to the disturb- 
ance of interest: a pressure disturbance occurring at a station just 

downstream of the steady-state position of the shock. This fixed sta- 
tion will be referred to as station 2. 


The relation between Ap-j_ and Apg can be determined frcm the con- 
tinuity and momentum equations of unsteady gas flow. For the quasi- one- 
dimensional flow of a nonviscous gas (with no external forces ), the equa- 
tions are 


Continuity: 


Momentum: 


ox 




du du 

dt + u d£ 



= 0 


Frcm these fundamental equations and the relation 


p = f r 

T T 
P Pr 


p„ (S-S_)/c, 


the following Bet of equations can be derived: 


| + (u + a) ^ 


It + (u - a) 


+ a 


- a 


du / , \ du 

dt + (u + a) SE 


C J 


du t \ du 

dt + (u - a) 


a 2 DS 
Dt 


a 2 DS 
c Dt 


a J DA 
A Dt 

af DA 
A Dt 


> 


(4) 
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Downstream of the moving shock, of course, each, particle of gas maintains 
a constant value of .its entropy (although different particles will have 
different entropies). The substantial derivative of entropy, therefore, 
vanishes : 


E6 _ 
Dt 


0 


In addition, If the area of the duct is not time dependent, the substan- 
tial derivative of area reduces as follows: 


DA SA 
Dt = U "Sic 


With these simplifications, equations (4) become 



(u + a) 




+ al 


5u 

St 


+ (u + a) 



+ (u - a) 


chi 

Su 


a 2 u 9A 

~ Sk 

a 2 u dA 
A Sx 


> 


(5) 


In linearizing equations (5), It is convenient to make an assumption 
that appears consistent with the quasi-one-dimens ional treatment. It is 
assumed that partial derivatives with respect to position x for the 
steady-state condition are of the order of the small-perturbation quanti- 
ties being considered. That is, the area variation of the duct is assumed 
to be gradual. Upon linearizing, therefore, for small perturbations from 
the steady-state condition, equations (5) became 


^ + (u - a) 


+ a 


- a 


, ouu. 

3T 


+ (,u 


a 2 u 5AA 
A Sx 

a 2 u dAA 
A dx 




( 6 ) 


where, once again, only the A quantities are time dependent. 


For the purposes of the present analysis, equations (6) may be 
greatly simplified. This simplification results from the fact that the 
bracketed quantities in these equations represent differentiation along 
the paths of the sonic waves. If the analysis, therefore, is limited to 
transients that are slow in relation to the transit times for sonic waves, 
the partial derivatives with respect to time, may be omitted from the wave- 
path derivatives. Such a limitation is not overly restrictive for the 
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problem considered herein. Therefore, with this limit imposed, equations 
(6) may be written as 

u + a 3An r \ SiAu a 2 u cSiAA. 1 

■ + a) = - — "sr 


p 

u - a 
P 


+ a( 


f \ &Au a^u dAA. 
^(u - a) 


(?) 


Multiplying the second of equations (?) by ^ and then adding the 

result to the first of equations (7) yields the following : 


2 (u + a) 
P 


a 2 u f ~ u + a\ &AA 

— r l, 1 + t^sJ sr 


(a) 


When equation (8) is multiplied by dx and no distinction is made between 
the steady-state conditions of stations 1 and 2, it can be integrated to 


ppap u; 

Api = ^>2 " “X 


"1 ut 


2 ^1 


( 9 ) 


which can be put into the following form: 

,2 


, A rPi m£ (£* o\ 

* ^2 - Ay M | . ! {tx 0 ) Be 0 


( 10 ) 


where the limit can now be defined, frequencywise, as 


CD « 


ap± u i 

(Ax^) 

^ u'max 


The relation between &p ± and Ap 2 having thus been determined, the 
equation relating shock position Ax 0 to the pressure disturbance of 
interest Ap g can now be formulated. This formulation is done by combin- 
ing equation (10) with equation (3); the manipulation yields 



TPi Ml f ^o A 

+ A o m| - iJWoy ss 


Ax 0 - 


4r d Oo> 

X + 1 a Q dt 


(H) 
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From, equation (ll) it can readily be verified that, in the Laplace 
domain, the first-order lag relation 


Ax 0 (s) 
APp(b) ~ k 


1 + as 


(12a) 


exists. The time constant a and the gain k are defined in terms of 
a dimensionless quantity 0 ’ as follows: 


a = 


v 0 


; °( 5 ) 

' u 'ss 


(12b) 


k = - 


r + i i l a o 


Ao 


4 r p 0 Mo % f^o\ 

V*oJ. 


ss 


(12c) 


where 


0 1 = 


^0 A o 

^ a 0 Aq 



r + i £i £o 

4 r po p 0 



(I2d) 


r + l p i Ao % 


*0 *0 - 1 


Equation (l2d) can be reduced to a simpler form by using the familiar 
functions of Mq that define the parameters __ A£/A q , P-j/pq, 


f &0 


Pq/Pq, and (ref. 3) and by considering the quantities | — j- 

A 7* 


and 



to be equal to the corresponding derivatives. The reduced form 


ss 

is as follows: 




1/2 




(l2e) 
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Thus, the quantity a' , which is called the dime nsionless shock t ime 
constant, is dependent only on the steady-state value of the shock Mach 
number Mq. A curve showing the variation of the dimen sionless time con- 
stant with shock Mach number is presented in figure 1. (in determining 
this curve, a value of 7/5 was used for y.) As can be seen from, the 
curve, c' varies only slightly for a wide range of Mach numbers; it is, 
therefore, a convenient parameter to use. 


With the curve of figure 1, it becomes a simple procedure to calcu- 
late the time constant a and the gain k that define the shock dynamics 
for a given duct and initial condition. As can be seen from equation 


(12b), the value of a will be positive for a diverging duct [ ~ - > 


and negative for a converging duct 


©4 


AA C 

These facts, of course. 


®4 


are consistent with observed stability phenomena of normal shocks and 
with the conclusions of reference 2. 


For the special case of a constant-area duct 



equations 


(12b) and (12c) show that both the shock time constant a and the gain 
k are infinite. For this case, therefore, a more useful form of equa- 
tion (12a) is as follows: 


Ax 0 (s) -L 

Ap 2 (s) “ 1 a 
k k 


which, for infinite a and k, gives 


Ax 0 (s) 

AP^Ts) = t 4r PqMq T 

~V + 1 a o J S 


In a constant-area duct, therefore, an integrator action exists between 
shock position and a downstream disturbance in pressure, the integrator 
rate being 


1 

4r 

r + 1 a 0 

In other words, the shock wave continues to move as long as the disturb- 
ance is present. This, again, is consistent with experimental 
observations . 
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CONCLUDING REMARKS 

According to the small -perturbation analysis, a downstream disturb- 
ance in pressure (within a specified frequency limit) causes the position 
of a normal shock in a duct to change in first-order lag manner. The 
time constant and gain of this lag are functions of steady- state condi- 
tions upstream of the shock, the configuration of the duct, and a dimen- 
sionless time constant that depends only on the steady-state value of the 
shock Mach number. For the special case of a constant-area duct, this 
lag reduces to an integrator action. 


Lewis Flight Propulsion Laboratory 

National Advisory Committee for Aeronautics 
Cleveland, Ohio, July 12, 1957 
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Shock Mach number, Mq 

Figure 1. - Variation of dimensionless shock time constant with shook Mach number. 
[Ratio of specific heats, 7/5.) 


i 


MCA TH 4090 


